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AN  ANALYTICAL  DESCRIPTION  OF  THE 
BEHAVIOR  OF  GRANULAR  MEDIA* 

by 

Jay  B.  Weidler**  and  Paul  R.  Paslay*** 

Normality  of  the  plastic  strain  increment  vector  to  the  loading  sur¬ 
face  is  not  always  observed  for  granular  media.  Beginning  from  a  thermo¬ 
dynamic  basis  and  considering  only  isothermal  processes,  a  series  of 
physically  appealing  assumptions  is  employed  to  increasingly  delimit  the 
possible  directions  of  the  plastic  strain  increment  vector.  The  results 
of  the  analysis  indicate  that  on  certain  portions  of  the  loading  surface 
normality  should  hold  while  on  the  remainder,  envelopes  or  fans  are  possi¬ 
ble.  Finally,  several  alternative  restrictions  are  offered  which  provide 
a  unique  stress-strain  relation  at  every  point  on  the  loading  surface. 
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Introduction 

A  granular  media  is  qualitatively  described  as  an  aggregate  of  particles 
whose  resistance  to  inelastic  flow  in  enhanced  by  an  increase  in  the  hydro¬ 
static  component  of  the  external  forces.  The  voids  of  such  a  material  are 
considered  to  be  interconnected  and  these  interstices  may  be  dry,  partially 
saturated  or  completely  filled  with  liquid.  It  will  be  presumed  for  the 
system  in  question  that  drainage  is  rapid  enough  so  that  the  pressure  in  the 
liquid  remains  constant  under  a  quasi-static  load  change.  The  particles  com¬ 
prising  the  aggregate  are  considered  as  rigid  bodies  of  unspecified  size  or 
shape.  The  over-all  size  of  the  particle  is  assumed,  however,  to  be  small 
enough  with  respect  to  the  size  of  the  system  to  permit  treatment  by  the 
methods  of  continuum  mechanics.  The  theory  to  be  presented  here  allows  for 
cohesion  between  particles,  and  it  presumes  time-independent  and  isothermal 
behavior. 

Historically  a  mathematical  formulation  of  the  division  between  attain¬ 
able  and  unattainable  stress  states  in  granular  media  was  proposed  first  by 
Coulomb  [1]*.  This  was  generalized  later  by  Mohr  [2].  The  form  of  the  line 
demarking  the  regions  is  shown  in  Figure  1  for  a  two-dimensional  (t  -  p) 
stress  space.  Until  recently  the  major  aim  of  those  working  in  the  mechanics 
of  granular  media  has  been  in  finding  states  of  stress  in  equilibrium  with 
the  applied  loads  for  generalized  inelastic  flow  conditions.  The  extent  of 
this  endeavor  is  described  in  a  recent  book  by  Sokolovskii  [3], 
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Numbers  in  brackets  designate  references  listed  in  the  bibliography. 
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In  contrast  to  the  above  mentioned  effort  the  literature  concerned  with 
the  generalization  of  the  stress-strain  relationship  had  to  be  deferred  until 
the  theory  of  plasticity  had  been  well-developed.  Generalizations  of  a  lin¬ 
earized,  two-dimensional,  Mohr-Coulomb  criterion  (shown  in  Figure  2)  to  in¬ 
clude  stress-plastic  strain  increment  relations  were  proposed  by  Drucker  and 
Prager  [4],  Shield  [5],  Drucker,  Gibson,  and  Henkel  [6],  Jenike  and  Shield 
[7],  and  Drucker  [8],  [9].  All  of  the  above  employed  the  concept  of  normality 
[10]  to  determine  the  direction  of  the  plastic  strain  increment  vector. 

Experimentally  it  is  observed  that  granular  media  are  inelastically  de¬ 
formed  almost  immediately  upon  application  of  stress  and  subsequently  work 
harden  until  the  stress  state  reaches  the  Mohr-Coulomb  envelope.  At  this 
point  if  the  material  was  initially  in  a  loose  state,  it  begins  to  flow,  while 
if  the  initial  state  was  dense,  work  softening  begins.  Although  much  effort 
has  been  devoted  to  experimentally  determining  the  shape  of  the  limiting  en¬ 
velope,  little  work  has  been  reported  describing  the  shape  of  the  intervening 
loading  surfaces.  Consequently  correlation  of  theory  to  experiment  has  rested 
largely  on  observation  of  the  direction  of  the  strain-increment  vector  at 
stress-states  on  the  limiting  surface.  In  Figure  3  if  normality  holds,  then 
the  strain  increment  vector  must  indicate  a  volume  increase.  Numerous  experi¬ 
ments  give  evidence  that  the  rate  of  volume  increase  lies  closer  to  the  direc¬ 
tion  of  zero  volume  change  than  to  the  direction  required  by  normality. 

Although  this  phenomenon  can  be  explained  through  the  use  of  corners  on 
the  loading  surface,  an  attempt  will  be  made  the  subsequent  discussion  to 
present  an  alternative  explanation.  In  granular  media  the  resistance  to  de¬ 
formation  is  furnished  by  two  types  of  friction.  The  first  of  these  arises 
on  the  contact  surfaces  of  adjoining  particles  termed  surface  friction.  The 
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second,  denoted  as  interlocking  friction  refers  the  resistance  offered  by 
the  interference  of  the  particles  themselves  to  changes  of  their  relative 
positions.  It  has  been  previously  pointed  out  that  normality  is  not  a  law 
of  nature  and  that  in  frictional  systems,  normality  does  not  hold  in  gen¬ 
eral  [11],  [12]. 

Perhaps  the  most  significant  alteration  incorporated  in  the  theory  pre¬ 
sented  here  is  the  dependence  of  the  loading  surface  on  the  specific  volume 
of  the  material.  This  phenomenon  has  been  discussed  by  Jenike  and  Shield 
[7].  If,  for  a  given  stress  state,  the  flow  makes  the  material  more  dense, 
then  work -hardening  will  generally  occur  while  if  the  flow  is  such  that  the 
volume  increases,  the  material  will  generally  work  soften.  Whenever  a  theory 
includes  this  feature,  the  direction  of  the  plastic  strain  increment  vector 
assumes  the  additional  significance  that  besides  determining  the  flow  state 
it  determines  the  subsequent  loading  surface. 

The  purpose  of  this  paper  then  is  to  develop  a  stress  —  plastic  strain 
increment  relation  which  accommodates  the  observed  phenomena. 

General  Considerations 


The  theory  to  be  presented  is  deduced  by  observing  the  restrictions  re¬ 
sulting  from  a  series  of  physically  motivated  assumptions.  These  assumptions 
are  presented  in  a  sequence  which  begins  from  a  thermodynamic  basis  and  con¬ 
tinues  to  the  point  where  a  definite  theory  can  be  established.  Each  of  these 
assumptions  implies  its  predecessors  but  at  the  same  time  the  plausibility  of 
each  assumption  decreases.  In  the  ensuing  development  tensile  stresses  are 
taken  as  positive. 
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The  theory,  as  was  previously  mentioned,  is  based  on  isothermal  behavior 
with  changes  of  state  accomplished  quasistatically.  The  state  of  the  material 
is  assumed  to  be  known  when  the  stress  and  plastic  strain  tensors  are  known. 
Elastic  strains  are  neglected  and  the  initial  state  is  one  where  the  stress 
and  plastic  strain  tensors  both  vanish.  A  loading  surface  is  assumed  to  exist 
that  encloses  all  stress  states  which  do  not  induce  plastic  flow. 

The  configuration  of  the  loading  surface  is  determined  by  the  state  and 
is  defined  by  a  relationship  of  the  form: 


(1) 


Here  the  are  also  the  components  of  the  total  as  well  as  the  plastic 

strain  tensors,  since  the  elastic  components  are  assumed  equal  to  zero.  The 
direction  and  sense  of  the  vector  in  nine  dimensional  strain  space  with  com¬ 
ponents,  de. .^,  are  uniquely  determined  by  the  stress  state  a,  .  It  is  pos- 
13  K± 

tulated  that  for  each  stress  state  on  the  loading  surface  it  is  possible  to 
induce  plastic  flow  and  for  each  de.^  vector  there  is  at  least  one  point  on 
the  loading  surface. 

Finally  all  functions ,  including  the  loading  surface ,  are  assumed  to 
possess  sufficient  continuity  for  the  ensuing  development.  Specifically  such 
notions  as  corners  on  the  loading  surface  are  excluded. 


Conjugate  Points  and  the  Clausius  Inequality 

The  behavior  of  any  material  must  be  such  that  the  laws  of  thermodynamics 
are  not  violated.  The  first  such  law  is  a  balance  equation  between  work,  heat 
flow,  and  internal  energy.  Since  the  measurement  of  the  heat  flow  in  a 
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granular  medium  is  exceedingly  difficult ,  the  first  law  does  not  play  a  cen¬ 
tral  role  in  the  development  of  an  analytical  theory  of  behavior.  The  first 
law  can  thus  be  regarded  as  an  auxiliary  equation  which  determines  the  heat 
flow  in  any  process.  As  the  state  of  the  system  is  defined  by  ck  ^  and 
and  the  internal  energy  is  a  state  variable ,  it  is  clear  that  the  internal 
energy  per  unit  current  volume,  U  ,  may  be  expressed  by 

U  =  U(o. .,£..)  (2) 

il  il 


In  general  one  would  anticipate  that  when  the  elastic  strain  is  neglected  U 
would  be  dependent  only  on  the  plastic  strain.  This  restriction  is  intro¬ 
duced  later.  Thus  the  first  law,  or  balance  of  energy  is  given  by 


dQ  =  dU  -  a.,  de..  (3) 

il  il 


where  dQ  is  the  differential  of  the  heat  flow  per  unit  current  volume  into 
the  system  and  de^_.  wor^  Per  unit  current  volume  done  on  the  sys¬ 

tem.  Since  de„  equals  de^^  (elastic  strains  being  neglected).  Equation  (3) 
may  be  rewritten  as 


dQ  =  dU  -  a. .  de. . 

il  il 


(4) 


This  result  gives  no  restrictions  on  the  behavior  of  granular  media.  It  only 
gives  a  means  for  finding  dQ  for  a  differential  change  of  state  in  an  element 
of  the  material. 

The  second  law  of  thermodynamics  may  be  conveniently  stated  in  terms  of 
a  cycle.  Il'yushin  [13]  has  used  this  concept  to  develop  a  thermodynamic  pos¬ 
tulate  for  elastic-plastic  materials.  A  cycle  is  defined  as  a  continuous 
change  of  state  where  the  initial  and  final  states  coincide.  Before  examining 
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the  consequences  of  the  second  law,  the  concept  of  a  cycle  in  a  granular  sys¬ 
tem  must  be  clarified.  The  requirement  of  returning  to  the  initial  state  in 
a  cycle  implies  a  zero  net  change  in  the  plastic  strain.  It  further  implies 
by  Equation  (1)  that  the  initial  and  final  loading  surfaces  are  the  same. 
Consequently  a  cycle  in  a  granular  system  requires  that  the  original  packing 
or  an  equivalent  packing  can  be  restored  or  established  by  a  properly  chosen 
sequence  of  stress  states.  Isotropically  work  hardening  materials  provide  an 
example  where  this  requirement  on  cycling  cannot  be  satisfied  by  isothermal 
processes .  Here  the  size  of  the  loading  surface ,  an  additional  state  vari¬ 
able,  always  increases  even  though  the  net  plastic  strain  can  be  made  equal 
to  zero.  In  reality,  permanent  changes  to  the  granular  system  do  occur  during 
plastic  deformation,  but  it  is  assumed  that  they  are  of  secondary  importance 
to  the  over-all  behavior.  Such  changes  include  modification  of  surface  fric¬ 
tion  factors,  breakage  or  abrasion  of  grains,  and  permanent  deformation  or 
rupture  of  cohesive  bonds. 

Having  discussed  the  concept  of  a  cycle  in  general  terms,  the  following 
thought  experiment  illustrating  the  concept  will  be  useful  in  the  subsequent 
development  of  the  theory. 

First  applying  a  loading 


o. . 

il 


(5) 


where  is  a  constant  stress  vector  and  k  is  a  parameter  which  has  a  zero 
value  in  the  reference  state  and  increases  monotonically  to  as  the  stress 

state  reaches  the  initial  loading  surface.  Hence 


o.©=  £>i.. 

ii  Lii 

F(o.  ©,  0)  =  0 
il 


(6) 


(7) 
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By  a  previous  postulate  the  direction  and  sense  of  the  strain  vector  are 
determined  by  ck  For  this  thought  experiment  a  differential  amount  of 
strain,  de„^,  is  induced  to  which  there  corresponds  the  nine  component 
vector,  di?3,  in  strain  space.  The  assumption  that  corresponding  to  this 
there  is  a  differential  change  in  stress  dc^  CD  is  substantiated  experimen¬ 
tally  for  most  soils  and  is  accepted  in  what  follows .  Consequently 


F(o.CD+  do.  CD,  de..^3)  =  0 

13  13 


defines  the  new  yield  surface.  In  order  to  complete  this  cycle  the  element 
must  be  returned  to  zero  stress  and  zero  strain  as  these  are  the  state  vari¬ 
ables.  Consequently  the  negative  of  the  existing  plastic  strain  must  be  in¬ 
duced.  A  previous  assumption,  insuring  the  availability  of  all  strain  vec¬ 
tors,  indicates  that  for  at  least  one  stress  state,  say  c^CD,  on  the  new 
loading  surface,  the  initial  change  of  plastic  strain  will  be  opposite  to 
de*3 .  Such  a  stress  state  is  said  to  be  the  conjugate  stress  state  to  CD. 

The  next  increment  of  the  loading  program  is  to  change  the  stress  state  from 

O.CD+  do.  CD  to  a.  CD  in  such  a  way  that  no  additional  plastic  flow  occurs. 

13  13  13 

It  is  assumed  that  such  a  loading  path  can  be  found.  At  this  point  a  differ¬ 
ential  change  of  strain  is  induced  which  eliminates  the  existing  plastic 
strain.  This  is  possible  since  the  initial  plastic  strain  increment  vector 
at  0^ CD  is  directed  opposite  to  d?3.  The  final  step  in  the  loading  program 
is  to  return  to  zero  stress  without  inducing  further  plastic  deformation. 

This  gives  the  same  initial  and  final  states  to  the  cycle;  i.e.,  =  e^.  =  0. 

This  cycle  is  presented  schematically  in  Figure  4.  The  above  cycle  of  load¬ 
ing  may  be  extended  to  include  finite  changes  of  strain  as  well  as  cycles  be¬ 


ginning  from  an  arbitrary  state, 
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It  is  now  of  interest  to  investigate  the  consequences  of  applying  the 
second  law  of  thermodynamics  to  such  cycles.  For  an  isothermal  process  the 
second  law  of  thermodynamics  requires  that 


dO  £  0 


(9) 


where  Vjy  implies  integration  around  a  closed  cycle  and  dQ  is,  as  mentioned 
before,  the  heat  flow  per  unit  current  volume  into  the  system.  Actually  what 
is  meant  by  an  isothermal  process  is  not  defined  since  temperature  has  not 
been  defined.  In  this  sense  it  may  be  more  appropriate  to  regard  Equation  (9) 
as  a  postulate.  Now  since  the  internal  energy  change  in  any  closed  cycle  is 
zero.  Equation  (9  )  and  Equation  (4)  yield 


ck..  de^13  %  0  (10) 

Thus  for  the  cycle  depicted  schematically  in  Figure  4,  Equation  (10)  may  be 
written  as 


a.®de..^  -  a..©de..^  >,  0 
il  ID  ID  ID 


(ID 


if  products  of  the  second  order  in  the  differentials  of  stress  and  strain  are 
neglected.  On  factoring  out  de..^,  the  result 


(a.®-  a.®)  de..^  >,  0 
ID  ID  ID 


(12) 


may  be  interpreted  as  requiring  the  dot  product  of  the  vector  (a^CD-  CSb 
and  the  vector  de.^  to  be  greater  than  or  equal  to  zero.  In  other  words  the 
projection  of  the  vector  de*5  on  (3S>  -  32>)  in  the  nine  dimensional  stress- 
strain  space  must  be  positive  or  zero.  A  schematic  representation  of  this  is 
shown  in  Figure  5  where  the  admissible  directions  of  de^  are  depicted. 
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If  for  t  two  conjugate  points,  ®  and  ck  £D  »  exist  the  limitations 

on  the  direction  of  d?  are  indicated  schematically  in  Figure  6.  Since  for 
each  stress  state  on  the  loading  surface ,  it  is  possible  to  induce  plastic 
flow,  the  conjugate  stress  states  must  be  located  so  that  a  plane  can  be 
passed  through  £D such  that  all  conjugate  points  are  in  the  plane  or  on  one 
side  of  it. 

A  Less  Restrictive  Concept  of  State 

It  is  clear  from  the  preceding  section  that  more  definite  assumptions 
must  be  made  before  the  direction  of  d?  can  be  uniquely  determined  for  a 
stress  point  on  the  loading  surface.  Now  any  sufficiently  restrictive  pos¬ 
tulate,  as  long  as  it  satisfies  Clausius'  Inequality,  Equation  (10),  can  es¬ 
tablish  a  definite  stress-strain  increment  relationship;  for  example,  see 
Drucker  [10].  The  approach  to  be  followed  in  this  paper  is  to  continue  to 
make  further  assumptions,  usually  in  thermodynamic  terminology,  finally 
reaching  a  point  where ,  for  each  point  on  the  loading  surface ,  de^  is 
uniquely  determined.  The  first  assumption  in  this  process  is  presented  in 
this  section. 

The  preceding  analysis  depended  on  a  closed  cycle  which  had  the  same 
initial  and  final  state .  In  the  case  of  granular  media  it  appears  to  be  of 
interest  to  explore  the  consequences  of  specifying  the  state  in  a  less  re¬ 
strictive  manner  than  was  done  in  the  preceding  section.  Let  the  state  be 
now  determined  if  and  v  ,  the  specific  volume,  are  known. 

Since  the  differential  change  in  v  is  given  by 

dv  =  v  de . .^ 

li 


(13) 
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and  as  the  configuration  of  the  loading  surface  must  be  the  same  for  equiva¬ 
lent  states,  i.e. 


F(o. . ,  v)  =  0 


(14) 


then  the  process  of  a  closed  cycle,  similar  to  that  of  Figure  4,  is  shown  as 
Figure  7.  Other  than  the  modified  definition  of  the  state  and  Equation  (14) 
replacing  Equation  (1),  the  preceding  theory  is  unaltered. 

In  this  case  it  is  anticipated  that  more  than  one  conjugate  point  exists 
for  every  point  on  the  loading  surface .  The  reader  should  note  that  at  a  con- 

i ©  p  .  p 


jugate  point  in  this  cycle  de^ 


-  de.  . 


0 


d£ii  ku"t  that,  in  general,  de„™  t 


In  fact  it  is  anticipated  that  any  point  on  certain  regions  of  the 


loading  surface  would  be  a  conjugate  point  for  cm.p. 

The  restriction  placed  on  the  closed  cycle  shown  in  Figure  7  by  Equation 
(10)  is 


1]  1] 
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(15) 


for  every  conjugate  stress  state 

It  is  clear  that  this  less  restrictive  concept  of  state  will  not,  by  it¬ 
self,  lead  to  a  definite  specification  of  de^.  Subsequent  assumptions,  how¬ 
ever,  will  incorporate  this  modification.  Owing  to  the  fact  that  the  results 
presented  in  this  section  are  a  special  case  of  the  preceding  section,  all  re¬ 
sults  are,  thus  far,  consistent. 


Discussion  of  the  Strain  Energy  Function 

In  the  previous  section  the  state  variables  were  reduced  to  and  v  . 
Implied  in  this  assumption  was  a  concept  that  the  only  permanent  deformation 
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altering  the  response  of  the  material  to  subsequent  changes  of  stress  was  a 
change  in  the  specific  volume. 

The  assumption  made  in  this  section  is  that  the  internal  energy,  U  ,  is 
a  function  of  the  specific  volume  only,  i.e., 

U  =  U(v)  (16) 

and 

dU  =  ^  dv  =  vde..P  (17) 

dv  dv  n 

Physically,  the  basis  for  this  assumption  lies  in  the  mechanisms  avail¬ 
able  in  granular  media  for  storing  internal  energy.  If  an  aggregate  of  par¬ 
ticles  decreases  in  volume,  it  is  possible  for  some  of  the  particles  to  be 
wedged  between  their  neighbors  in  such  a  manner  that  on  a  release  of  stress, 
residual  forces  remain  in  the  system.  Should  subsequent  stressing  alter  the 
existing  volume  then  it  is  plausible  that  these  residual  forces  could  be  aug¬ 
mented  or  released. 

To  place  some  bounds  on  the  possible  variations  of  U  with  v  ,  the 
pressure-volume  curve  typical  to  granular  media  is  shown  as  Figure  8.  This 
curve  provides  a  limit  on  — ,  i.e. 

0  5  -  —  v  $  p*  (18) 

dv 

so  that  Figure  9  results,  p*  represents  that  value  of  p  where  the  stress 

path  -  p6..  intersects  the  current  loading  surface.  Now  since  p*,  the  hydro- 
J 

static  pressure  is  positive,  ^  must  be  less  than  or  equal  to  zero  and  hence 
the  internal  energy  —  specific  volume  curve  has  a  shape  as  depicted  in  Fig¬ 


ure  10. 
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With  this  assumption  on  the  internal  energy,  it  is  of  interest  to  apply 
a  load  to  the  system  causing  a  differential  change  of  strain  de.j^P  and  con¬ 
sequently  a  change  in  specific  volume  v  .  As  ^  is  anticipated  to  be  nega¬ 
tive,  there  will  be,  by  Equation  (17),  energy  released  from  storage  if  dv 
is  positive  and  energy  stored  if  dv  is  negative.  During  this  process  the 
consequences  of  requiring  that  the  sum  of  the  work  done  by  the  stress  field 
and  the  energy  released  from  storage  be  positive  semi-definite  may  be  ex¬ 
pressed  as 

dei;.P  -  dU  *  0  (19) 


or  using  Equation  (18)  as 


(oij 


dU 


v6. . )  de. .  i  0 


dv  lj 


13 


(20) 


A  schematic  representation  of  Inequality  (20)  is  shown  as  Figure  11.  In 
this  figure  the  reader  should  note  that  p  =  p*  serves  to  govern  the  magnitude 
of  the  vector  “  v  as  0  £  v  ^  $  p*.  This  requirement  is  a  result  of  a 
postulate  concerning  the  way  in  which  energy  is  stored  and  released  in  a  soil 
and  has  no  deeper  meaning  than  this.  Comparison  of  Equation  (4)  and  Inequal¬ 
ity  (20)  shows  the  equivalence  of  this  requirement  to  that  obtained  by  postu¬ 
lating  that  during  isothermal  flow  the  heat  flow  from  the  soil  is  positive  — 
semi-definite. 

If  Inequality  (20)  is  written  for  cr^CD,  where  is  a  stress  state  on 
the  loading  surface,  then  for  a  differential  change  in  plastic  strain  de^^ 
(recall  v  >  0) 


13  dv  11 


(21) 
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Now  for  a.  ©a  stress  state  which  is  a  conjugate  to  a.®.  Inequality  (20) 

/T\  ^ 


can  be  written  for  a  change  in  plastic  strain  de^ .  such  that 


dc.A  -  de..-® 
11  11 


(22) 


as 


il  il 


*  0 


(23) 


Inequalities  (21)  and  (23)  imply 

cr..®de..^"^  +  a.@de..^^^  0  (24) 

il  il  il  il 

which  is  equivalent  to  Inequality  (15).  Consequently  the  restriction  of  In¬ 
equality  (20)  implies  all  previously  derived  restrictions. 

As  the  direction  for  d?3  is  still  not  uniquely  determined  the  following 
section  introduces  further  limitations  to  be  imposed  on  the  soil  behavior. 


Stability  Considerations 

A  direct  postulate  for  stability  of  soil  would  be  unreasonable  as  experi¬ 
ments  indicate,  that  soil  is,  in  the  usual  sense,  unstable  for  certain  kinds 
of  loading.  However,  it  seems  hopeful  to  try  to  make  some  limitation  on  the 
material  behavior  corresponding  to  a  stability  postulate.  As  previous  reason¬ 
ing  has  been  based  on  a  storage  of  energy  (change  of  internal  energy)  deter¬ 
mined  by  the  change  of  specific  volume ,  it  is  appealing  to  introduce  this  con¬ 
cept  into  a  "stability"  postulate. 

One  approach  which  leads  to  a  definite  restriction  is  to  postulate  some¬ 
thing  about  the  energy  balance  when  a  set  of  loads  is  superimposed  on  an 
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initial  set  of  loads.  Presume  that  the  initial  loading  induces  the  stress 
state  a^j*«  Obviously  this  initial  stress  must  lie  in  or  on  the  loading 
surface  corresponding  to  the  initial  state .  Now  the  superimposed  stress 
state  is  applied  so  that  the  total  stress  is  at  some  point  on  the  yield 
surface.  If  a  differential  change  of  plastic  strain,  de^..^,  is  induced  then 
it  is  postulated  that  the  sum  of  the  work  done  by  the  superimposed  stresses , 
a..  -  a..*,  during  the  change  and  a  fraction,  a,  of  the  energy  released  be 
greater  than  zero.  When  energy  is  stored  due  to  a  volume  decrease,  or  when 
the  volume  is  unaltered,  it  is  only  required  that  the  work  done  by  the  arbi¬ 
trary  superimposed  stresses  be  positive-semi-definite.  This  leads  to 


(a..  -  a..*)de..P  +  a  ^  v  de.  z  0 
i]  i]  \  dv  n/ 


iD  il  il 

where  the  notation 


<*>  ■* 

if 

O 

A 

X 

<*>  - 0 

if 

x  s  0 

is  used  and 


0  $  a  $  1 


As  the  point  maY  be  any  point  on  or  inside  the  yield  surface,  -  c„*) 

may  be  interpreted  as  a  vector  whose  tail  is  in  or  on  the  yield  surface  and 
whose  tip  is  on  the  yield  surface.  Anticipating  that  ^  is  negative  the  sta¬ 
bility  requirement  may  be  rewritten  in  the  following  form 

.  .  .  dU  X£kk/^J  p  . 

(oij  ~  °ij5>  -  6ij  ad7  v  — T)deij  *  0  (28) 

dekk 


which  yields  the  graphical  interpretation  shown  in  Figure  12. 
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When  0  the  stability  requirement  Inequality  (28)  simply  requires 

the  projection  of  de..p  on  ( cr . .  -  a..*  -a  4^- v  6..)  to  be  positive. 

F  J  13  i3A  13  dv  13 

When  de..P  <  0  the  above  requirement  reduces  to  (a..  -  a..*)  de-.P  >,  0. 

ii  i3A  13  13 

In  other  words ,  referring  to  Figure  12 ,  the  boundary  of  the  shaded  region 
coincides  with  the  loading  surface.  Consequently  when  de.j^P  <  0  the  stability 
requirement  leads  to  normality  of  d?5  with  the  loading  surface.  The  sense  of 
de^  is  outward. 

Further  when  de..P  <  0  at  point  P  on  the  limit  surface  then  (a..  -  a..*). 


11 


i3p  13 


de„P  £  0  implies,  in  geometrical  terminology,  that  all  points  enclosed  by  the 
loading  surface  are  on  one  side  of  the  tangent  plane  to  the  loading  surface  P. 
This  latter  property  will  be  referred  to  as  convexity. 

Now  consider  point  A  in  Figure  12.  Either  de.^  <  0  or  not.  Assume  that 
d&  J5  <  0,  then  outward  normality  must  hold.  According  to  Figure  12  if  de*5 
were  normal  to  the  loading  surface,  it  would  have  a  projection  on  the  = 
-p<5.^  line  in  the  direction  of  increasing  tensile  stress  so  that  de^.^  i  0 


contrary  to  the  assumption.  Consequently  at  point  A  d£^..P  £  0.  The  admissi¬ 
ble  directions  for  de^  deduced  from  Inequality  (28)  are  shown  in  Figure  13. 

Point  B  in  Figure  12  is  inside  the  shaded  area  if  de..^P  i  0.  Clearly  no 
non-trivial  deP  is  possible  such  that  Inequality  (28)  is  satisfied  if  de..^  %  0 
so  that  de^P  $  0.  It  had  previously  been  assumed  that  plastic  flow  could  be 
induced  for  every  point  on  the  loading  surface.  Consequently  outward  normality 
of  deP  to  such  points  as  B  must  exist. 

Finally  consider  point  C.  If  de^p  <  0  is  assumed  the  outward  normality 
condition  of  deP  to  the  loading  surface  results  whereas  if  de.^  i  0  is  assumed 
the  conditions  on  the  direction  of  dep  deduced  from  Inequality  (28)  and  shown 
in  Figure  14  result.  Consequently  either  de^P  i  0  or  de..P  <  0  may  occur  at 


ii 


point  C. 
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It  is  now  of  interest  to  determine  under  what  conditions  Inequality  (28) 
implies  satisfaction  of  Inequality  (20).  In  the  case  that  the  volume  in¬ 
creases  the  two  inequalities  are  equivalent  if 


9U  ( 1  \  r 

v  37  u-a)hl 


(29) 


If  the  volume  decreases  the  two  inequalities  are  equivalent  if 


8U 

V  -r —  0.  . 
3v 


(30) 


In  other  words,  satisfaction  of  Inequality  (28)  implies  satisfaction  of  In¬ 
equality  (20)  if  the  two  values  for  given  by  Equations  (29)  and  (30)  are 

3U 

admissible  (i.e.  inside  the  loading  surface).  Since  0  £  -v  -r—  $  p*  the  re- 

oV 

strictions  do  not  seem  unreasonable  and  consequently  will  be  assumed  satisfied 
in  the  following.  This  means  then  that  satisfaction  of  Inequality  (28)  auto¬ 
matically  implies  satisfaction  of  all  previous  requirements  which  have  been 
presented,  i.e..  Inequalities  (10),  (15)  and  (20) 

The  above  reasoning,  thermodynamic  and  otherwise,  imposes  restrictions  on 
the  kind  of  incremental  stress-strain  relationship  that  can  be  used  with  a 
specified  loading  surface.  In  the  case  where  a  is  identically  zero  the  con¬ 
ditions,  similar  to  plasticity,  show  that  normality  between  d?  and  the  corre¬ 
sponding  tangent  plane  to  the  loading  surface  must  be  specified  everywhere. 

This  result  is  significant  because  stipulation  of  the  loading  surface  then 
automatically  gives  the  incremental  stress-strain  relation.  In  the  present 
case  where  a  is  less  than  zero  the  incremental  stress-strain  relation  is 
determined  for  some  points  on  the  loading  surface  (point  B  types)  and  delimited 


for  the  remainder. 
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The  result  of  the  above  restrictions  then,  unlike  plasticity,  is  that 
further  restraints  must  be  introduced  before  a  definite  direction  for  d?3 
corresponding  to  each  stress  point  on  the  loading  surface  can  be  established. 
In  order  to  be  of  practical  use  in  solving  problems  a  definite  choice  of  di¬ 
rection  for  d?3  at  each  point  on  the  loading  surface  must  be  established. 

In  attempting  to  introduce  further  postulates  concerning  soil  behavior 
those  postulates  compatible  with  the  previous  set  of  assumptions  should  be 
considered  first.  That  is,  it  would  be  appealing  to  insure  that  the  previous 
set  of  postulates  could  be  satisfied.  The  next  section  presents  several  pos¬ 
sible  restrictions  which  lead  to  definite  choices  for  d?3  which  are  compatible 
with  the  theory  thus  far  developed. 

Possible  Restrictions  Leading  to  a  Definitive  Incremental  Stress-Strain 
Relationship 

The  three  restrictions  to  be  presented  in  this  section  which  lead  to  de¬ 
finitive  incremental  stress-strain  relationships  are  so  chosen  that  the  pre¬ 
ceding  theory  is  automatically  satisfied.  This  is  accomplished  by  only  con¬ 
sidering  directions  of  de^  admissible  from  the  preceding  theory  (hereafter 
called  admissible  directions)  and  then  making  a  statement  regarding  them  which 
leads  to  a  unique  choice  of  d?  at  each  point  on  the  loading  surface.  It  is 
convenient  to  define  the  length  of  a  d?3  vector  to  be 

length  of  d?3  =  jde^j  =  ^de^^  de^*3  (31) 

The  first  choice  presented  can  be  conveniently  stated  in  the  following 
proposition : 
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"Of  all  the  admissible  directions  of  d?  at  a  point  on  the  loading  sur¬ 
face  the  actual  de^  maximizes  the  plastic  work  a^de^^,  done  by  the  external 
loading  when  compared  to  all  other  admissible  de^  of  equal  length."  Figure  15 
shows  schematically,  the  consequences  of  such  a  proposition  (d?3  vectors  marked 
by(l)).  The  points,  B,C,A,A’,A"  correspond  to  points  B,C,A,A,A  of  Figures  12, 
13,  and  14.  This  choice  of  de*3  was  found  by  determining  which  projection  of 
de*3  on  the  vector  o  was  a  maximum  and  yet  lay  within  the  fans  provided  by  the 
previous  section. 

The  second  alternative  is  expressed  by  the  postulate 

"Of  all  the  admissible  de*3  at  a  point  on  the  loading  surface,  the  actual 

o . .de . .P 

de*3  maximizes  the  expression  {— ^ — — —  }  when  compared  to  all  other  admissible 

de, ,  ^ 

Kk 

de  . "  The  quantity  maximized  here  is  seen  to  be  the  ratio  of  the  plastic  work 
and  the  volume  change.  The  de*3  determined  by  this  postulate  are  shown  again  in 
Figure  15  by  the  vectors  marked  (2). 

The  final  postulate  presented  is  expressed  by  the  proposition: 

"Of  all  the  admissible  vectors  d'e*3  at  a  point  of  the  loading  surface  the 
actual  de*3  maximizes  the  dissipation  (o.^  -  ^  v  6^  )de„^  when  compared  to  all 
other  admissible  de*3  of  equal  length."* 

This  result  is  shown  again  on  Figure  15  by  the  vectors  marked© 

Which  of  these  three  postulates  and  possibly  some  others  not  presented  is 
the  correct  one  can  only  be  answered  by  further  analytical  and  experimental 
work.  The  validity  of  all  the  previous  assumptions  as  to  their  applicability 
to  problems  of  granular  media  can  again  only  be  answered  by  further  effort. 


This  formulation  was  pointed  out  to  the  authors  by  Professor  D.  C.  Drucker 
in  a  private  communication. 
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Discussion : 

In  order  to  improve  upon  the  currently  used  plastic  deformation  theories 
of  granular  media,  a  theory  is  presented  here  which  attempts  to  overcome  some 
of  the  known  deficiencies  in  previous  theories.  Clearly  the  development  of  a 
more  sophisticated  theory  almost  always  brings  mathematical  difficulites  so 
that  a  balance  must  be  maintained  between  sophistication  and  tractability. 

The  authors  feel  that  adoption  of  an  appropriate  specific  loading  surface 
function  in  the  framework  of  the  above  theory  will  lead  to  the  ability  to 
solve  several  technically  important  problems  [14],  [15]. 

Several  assumptions  made  in  the  theory  should  be  verified  experimentally 
by  comparison  of  predicted  and  measured  results  for  technically  important 
problems.  Direct  verification  of  individual  assumptions,  such  as  the  sta¬ 
bility  postulate  used  here ,  is  not  usually  possible .  Often  confidence  in 
such  postulates  can  be  strengthened  by  finding  relatively  simple  mechanical 
models  for  the  material  which  satisfy  the  postulate.  Although  no  such  model 
has  yet  been  conceived ,  except  for  a  =  0 ,  the  authors  feel  that  the  postulate 
may  nevertheless  have  meaning  when  0  <  a  $  1. 


-21- 


Bibliography 


1.  C.  A.  Coulomb,  "Essais  sur  une  application  des  regies  des  maximis  et 

minimis  a  quelque  problems  de  statique  relatifs  a  1* archi¬ 
tecture,"  Mem.  Acad.  Roy.  Pres.  Divers  Savants,  Vol.  7, 

Paris,  1776. 

2.  0.  Mohr,  Abhandlungen  aus  dem  Gebiete  der  Technischen  Meckanik,  Wilhelm 

Ernst  und  Sohn,  Berlin,  1905. 

3.  V.  V.  Sokolovskii,  Statics  of  Granular  Media,  Pergamon  Press,  Inc., 

New  York,  1965,  Translation  by  J.  K.  Lusher,  English  transla¬ 
tion  edited  by  A.  W.  T.  Daniel. 

4.  D.  C.  Drucker  and  W.  Prager,  "Soil  Mechanics  and  Plastic  Analysis  on 

Limit  Design,"  Quarterly  of  Applied  Mathematics,  Vol.  10,  1952, 
pp.  157-165. 

5.  R.  T.  Shield,  "On  Coulomb's  Law  of  Failure  in  Soils,"  Journal  of  the 

Mechanics  and  Physics  of  Solids,  Vol.  4,  1955,  pp.  10-16. 

6.  D.  C.  Drucker,  D.  J.  Henkel,  and  R.  E.  Gibson,  "Soil  Mechanics  and  Work- 

Hardening  Theories  of  Plasticity,"  ASCE  Transactions,  Vol.  122, 
1957,  pp.  338-346. 

7.  A.  W.  Jenike  and  R.  T.  Shield,  "On  the  Plastic  Flow  of  Coulomb  Solids 

Beyond  Original  Failure,"  Journal  of  Applied  Mechanics,  Vol.  26, 
1959,  pp.  599-602.  '  ‘  ~  ' 

8.  D.  C.  Drucker,  "On  Stress-Strain  Relations  for  Soils  and  Load  Carrying 

Capacity,"  Proceedings  of  the  First  International  Conference 
on  the  Mechanics  of  Soil-Vehicle  Systems,  Turin,  June  1961, 
Edizioni  Minerva  Tecnica,  pp.  15-23. 

9.  D.  C.  Drucker,  "Concepts  of  Path  Independence  and  Material  Stability  for 

Soils,"  Proceedings  of  the  International  Symposium  of  the  Inter¬ 
national  Union  of  Theoretical  and  Applied  Mechanics,  Grenoble, 
France,  April  1-8,  1964. 

10.  D.  C.  Drucker,  "A  More  Fundamental  Approach  to  Plastic  Stress-Strain  Re¬ 

lations,"  Proceedings  of  the  First  U.  S.  National  Congress  of 
Applied  Mechanics,  ASME ,  1951,  pp.  487-491. 

11.  D.  C.  Drucker,  "Coulomb  Friction,  Plasticity  and  Limit  Loads,"  Journal  of 

Applied  Mechanics,  Vol.  21,  1954,  pp.  71-74. 

12.  D.  C.  Drucker,  "Some  Implications  of  Work  Hardening  and  Ideal  Plasticity," 

Quarterly  of  Applied  Math.,  Vol.  7,  No.  4,  (1950),  pp.  411-418. 


XVLT;  :■ .  S13 

ABERDEEN  EJOTOfc  SRGtJSD,  W& 


c-rp  -p  r  rj  rm* 


-22- 


13.  A.  A.  Il'yushin,  "On  the  Postulate  of  Plasticity,"  Prikladnaya  Matematics 

i  Mekhaniko,  Vol.  25,  1961,  pp.  503-507. 

14.  P.  R.  Paslay  and  Jay  B.  Weidler,  "The  Triaxial  Test  in  Cohesionless  Soil 

Mechanics,"  to  appear  as  a  Brown  University  Report. 

15.  Jay  B.  Weidler  and  P.  R.  Paslay,  "The  Effect  of  an  Elastic  Membrane  on 

the  Stability  of  a  Triaxial  Test  of  a  Time  Independent  Soil," 
to  appear  as  a  Brown  University  Report. 


mlw 


FIGURE  CAPTIONS 


Figure  1  - 


Figure  2 


Figure  3  - 

Figure  4 


Figure  5 


Figure  6 
Figure  7  - 

Figure  8 


Figure  9 


Figure  10  - 
Figure  11  - 


Division  Between  Attainable  and  Unattainable  Stress  States  in 
a  Granular  Media 

Two-dimensional  Linearized  Form  of  the  Mohr-Coulomb  Yield 
Criterion 

Observed  Discrepancy  in  Direction  of  d?^  on  Limit  Surface 

Schematic  Diagram  Illustrating  Concept  of  a  Loading  Cycle  in 
Stress-Strain  Space 

Schematic  Diagram  Showing  How  ^o^de^^  ^  ®  Restricts  the 
Directions  of  d?  in  Terms  of  the  Stress  o.  and  One  of 
its  Conjugate 

Case  of  Two  Conjugate  Stress  States 

Schematic  Diagram  Showing  Loading  Cycle  in  Stress-Strain  Space 
for  Material  Whose  State  is  Given  by  o. . ,  v 

Typical  Specific  Volume  Versus  Hydrostatic  Pressure  Curve  for 
Granular  Media 

Bounds  on  ^  v  Provided  by  Hydrostatic  Pressure  —  Specific 
Volume  Curve 

Dependence  of  Internal  Energy  on  Specific  Volume  if  ^  $  0 
Admissible  d?  Which  Satisfy  Inequality  (21) 


Figure  12 
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Figure  14  -  Determination  of  Admissible  Directions  for  d?*  at  Point  C. 

In  this  case  de^^  $  0  Leads  to  Outward  Normality.  Either 
de„^  *  ®  or  ^eij^  >  ®  May  Satisfy  the  Condition  Imposed 
by  Inequality  (29)  (for  ^  <  0 ) 

Figure  15  -  Three  Possible  Restrictions  Leading  to  Definite  Incremental 

Stress-Strain  Relationships.  (1)  Maximum  Work,  (2)  Maximum 
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FIG.  I  DIVISION  BETWEEN  ATTAINABLE  AND 

UNATTAINABLE  STRESS  STATES  IN  A  GRANULAR 
MEDIA 


FIG.  2  TWO-DIMENSIONAL  LINEARIZED  FORM  OF  THE 
MOHR-COULOMB  YIELD  CRITERION 


FIG.  3  OBSERVED  DISCREPANCY  IN  DIRECTION  OF  d€P 
ON  LIMIT  SURFACE 


FIG.  4  SCHEMATIC  DIAGRAM  ILLUSTRATING  CONCEPT  OF 
A  LOADING  CYCLE  IN  STRESS  -  STRAIN  SPACE 
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FIG.  5  SCHEMATIC  DIAGRAM  SHOWING  HOW 

<^crjjd€jjP>0  RESTRICTS  THE  DIRECTIONS  OF 

d?P  IN  TERMS  OF  THE  STRESS  cr©  AND  ONE 
OF  ITS  CONJUGATES  cry® 


FIG.  7  SCHEMATIC  DIAGRAM  SHOWING  LOADING  CYCLE 
IN  STRESS-STRAIN  SPACE  FOR  MATERIAL 
WHOSE  STATE  IS  GIVEN  BY  crjj  ,  v 
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FIG.  8  TYPICAL  SPECIFIC  VOLUME  VERSUS  HYDROSTATIC 
PRESSURE  CURVE  FOR  GRANULAR  MEDIA 


SPECIFIC 

VOLUME 


FIG.  9 


v 

SPECIFIC 

VOLUME 
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BOUNDS  ON  ~v  PROVIDED  BY 
d  v 

HYDROSTATIC  PRESSURE -- SPECIFIC 
VOLUME  CURVE 


DEPENDENCE  OF  INTERNAL  ENERGY  ON 
SPECIFIC  VOLUME  IF  4^  <  0 


FIG.  II  ADMISSIBLE  d€P  WHICH  SATISFY  INEQUALITY  (21) 


LOADING  SURFACE  SHIFTED  BY 


FIG.  12  A  SCHEMATIC  WAY  OF  DETERMINING  ADMISSIBLE 

(o-::  -CTii*-av4^  Bn)  VECTORS.  ANY  VECTOR 
ija  U  dv  'J 

ORIGINATING  IN  THE  SHADED  REGION  IS 

ADMISSIBLE  .  (  FOR  -7-  <  0) 
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FIG.  13  DETERMINATION  OF  ADMISSIBLE  DIRECTIONS  FOR  d€P 

AT  POINT  A.  (  FOR  4^  <  °> 

d  v 


FIG.  14  DETERMINATION  OF  ADMISSIBLE  DIRECTIONS  FOR  dep 
AT  POINT  C.  IN  THIS  CASE  denp<0  LEADS  TO 
OUTWARD  NORMALITY.  EITHER  d€jjp<0  OR  d€jjp  >  0 
MAY  SATISFY  THE  CONDITION  IMPOSED  BY 
INEQUALITY  (28)  (FOR  ^-  <  0) 
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FIG.  15  THREE  POSSIBLE  RESTRICTIONS  LEADING  TO 
DEFINITE  INCREMENTAL  STRESS  -  STRAIN 
RELATIONSHIPS.  (I)  MAXIMUM  WORK ,  (2)  MAXIMUM 
OF  WORK  DIVIDED  BY  VOLUME  CHANGE, 

(3)  MAXIMUM  DISSIPATION  (FOr4^<0) 
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